Phase Transitions of Charged Scalars 
at Finite Temperature and Chemical Potential 



Rachel A. Rosen 

Oskar Klein Centre for Cosmoparticle Physics, Department of Physics 
Stockholm University, AlbaNova SE-10691, Stockholm, Sweden 

Abstract 

We calculate the grand canonical partition function at the one-loop level for scalar 
quantum electrodynamics at finite temperature and chemical potential. A classical 
background charge density with a charge opposite that of the scalars ensures the 
neutrality of the system. For low density systems we find evidence of a first order 
phase transition. We find upper and lower bounds on the transition temperature 
below which the charged scalars form a condensate. A first order phase transition 
may have consequences for helium-core white dwarf stars in which it has been argued 
that such a condensate of charged helium-4 nuclei could exist. 



1 Introduction and Summary 



There are two ways in which bosons can condense: via Bose-Einstein condensation or via 
spontaneous symmetry breaking. Bose-Einstein condensation occurs due to a conserved 
current associated with the bosons. Since the total number of particles (or other conserved 
quantity) can then be fixed, at low temperatures and high densities energy levels become 
"overcrowded." As a result the charge of the system must be stored the zero-momentum 
ground state. In theories with spontaneous symmetry breaking, a conserved current is 
not responsible for the condensation. Instead, the interactions of the bosons make it 
energetically favorable for them to condense. 

In this paper we consider the first mechanism, Bose-Einstein condensation, in the case 
that there is both a global U (1) symmetry corresponding to a conserved scalar number 
density as well as a local gauge symmetry. In particular, we consider scalar quantum 
electrodynamics in which the conserved scalar current allows us to introduce a nonzero 
chemical potential for the scalars. To ensure the neutrality of the system we consider 
scalars that reside in a constant background charge density of the opposite charge. In this 
work we treat this density as a classical, external background. 

To study the phase transitions of the system, we calculate the grand canonical partition 
function at the one-loop level. We find the critical temperature T c \ below which the ground 
state of the system is macroscopically occupied. Below this temperature some fraction of 
the scalar charge must reside in a condensate. We find also an implicit expression for the 
critical temperature T C 2 above which all of the scalar charge is stored in thermally excited 
states. In between these two temperatures both condensed and uncondensed solutions 
exist. 

For the usual neutral Bose-Einstein condensate these two temperatures coincide. At 
the critical temperature T c the expectation value of the scalars v which is the order pa- 
rameter of the phase transition goes continuously from zero to nonzero as the temperature 
of the system is lowered. This is characteristic of a second order phase transition. For the 
charged condensate we find markedly different behavior depending on the relative values 
of the number density n and the mass m. For low densities, when n 1 / 3 < a cm m, we find 
evidence of a first order phase transition: in between the critical temperatures T c \ and T C 2 
the value of v jumps discontinuously from zero. The lower critical temperature T c \ coin- 
cides with the critical temperature T c for the neutral condensate, while the upper critical 
temperature T C 2 can be several times larger than T c in the low density regime. This implies 
that for the charged condensate, condensation can occur at higher temperatures than for 
the neutral condensate. For high densities, when n 1 / 3 3> « em m, then T C 2 approaches T c \ 
and the transition appears to be quite nearly second order. 

The latter behavior can be understood through simple arguments. When the average 
kinetic energy of a boson is much greater than its Coloumb energy, the gas of bosons can 
be treated as nearly ideal. For this to be the case near the transition temperature, it is 
necessary that n 1 / 3 S> a cm m. Thus the physics of the charged condensate in the high 
density regime should be well approximated by the ideal, neutral condensate. The main 
interest of this paper is the low density regime, when n 1 / 3 < a cm m. In this regime the 
mass of the gauge boson in the condensate m 7 ~ et) becomes comparable with or greater 
than the transition temperature T c and a first order phase transition is possible. 

Such a system of charged scalars has been considered recently in the context of helium- 
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core white dwarf stars [H El El El [5] • In these works helium-4 nuclei play the role of the 
scalars while a degenerate electron gas is the background charge density. The number 
density of the nuclei n ~ (0.1 MeV) 3 is such that the average interparticle separation falls 
between the nuclear and atomic scales. Thus the nuclei are unable to form neutral atoms, 
while nuclear effects are negligible. It is also the case that the number density satisfies 
ra 1//3 < a em m, as m — 4 GeV. While standard theory holds that white dwarfs with carbon 
or oxygen cores will crystalize as they cool, it was argued in Refs. [21 [5] that for helium- 
core white dwarfs quantum effects will become significant before the crystallization point 
is reached. Thus instead of crystallizing, the helium nuclei can condense. 

In these works the condensation of helium-4 nuclei was described using a low energy 
effective field theory. The work of this paper is relevant to the helium-core white dwarfs 
if we take scalar QED to be a relativistic effective field theory. The relativistic theory is 
not the most appropriate description of the helium-4 condensate. It is overly restrictive in 
that it enforces Lorentz invariance - a symmetry we do not expect the low energy system 
of helium-4 nuclei and electrons to preserve. In addition it contains a heavy mode which 
we would expect to be beyond the scope of a low energy theory. However, it captures 
the significant features of the condensate including equivalent dispersion relations for the 
relevant degrees of freedom. 

It was recently argued in Ref. [6] that, in a neutral system of deuterons and electrons, at 
sufficiently high densities and low temperatures the spin-1 deuterons will condense. Such a 
system may be relevant to certain low-mass brown dwarf stars in which it is expected that 
such matter can exist in localized regions. In addition, because deuterons can condense 
at lower densities than helium, it may be possible to produce such high-density deuteron 
matter in laboratories in future shock wave compression experiments. The methods of 
this paper should be applicable to studying the phase transitions of these systems as well. 

Let us make a few comments on the literature. Gauge theories at finite temperature and 
symmetry breaking in gauge theories at finite temperature have been studied extensively, 
early works being El [9] (see also [TO] and references therein). In Ref. [11] spontaneous 
symmetry breaking in the abelian Higgs model at finite temperature was studied. It 
was found that a first order phase transition could take place if A <C g 4 where A is the 
coupling of the ^ 4 interaction and g is the charge of the scalar. For A > g 4 the phase 
transition is second order. The effect of a nonzero background fermionic current density 
at zero temperature was also considered in this work. It was found that an increase of the 
background charge density could lead to an increase in symmetry breaking. 

In Ref. [12] the abelian Higgs model with a background fermionic current density 
at finite temperature was studied in the limit that A > j 4 . The same model was also 
considered in Ref. [13] using a somewhat different method: instead of introducing a 
background fermionic current, a chemical potential associated with the conserved scalar 
current was introduced. The A > limit was taken here as well. Taking this limit allows 
one to ignore terms that mix the scalar and gauge bosonic fields. In our work we assume 
that the </> 4 interaction is negligible compared to the electrostatic interactions. Thus this 
mixing term plays a significant role, changing the spectrum of the theory in the condensed 
phase. Moreover, the emphasis of these two works is the effect of a conserved current 
on spontaneous symmetry breaking, rather than the Bose-Einstein condensation of the 
scalars due to the conserved current. The latter is the focus of this work. 
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The effects of nonzero chemical potentials associated with conserved charges on sym- 
metry breaking at finite temperature were also considered in Ref. [H]. In Ref. [15] 
an exact expression for the one-loop effective potential for cf> theory with a global U(l) 
charge at finite temperature was obtained. Refs. |16|. [IT] consider the screening of electric 
charge and photon polarization in scalar QED at finite temperature and density. Phase 
transitions in scalar QED with a nonzero chemical potential at zero temperature due to 
external magnetic fields were considered in [3]. 

The organization of this work is as follows. In the following section we briefly review 
the condensation of charged scalars at zero temperature. A more thorough discussion can 
be found in the recent review [18]. In section 3 we compute the grand canonical partition 
function at the one-loop level. We also treat the UV divergences of the theory. In section 
4 we find the critical temperature T c \ below which the scalars have condensed and we find 
an implicit expression for T C 2 above which all of the scalars are in thermally excited states. 
We plot the expectation value of the scalars v as a function of temperature for various 
values of the mass m and number density n of the scalars. We discuss the evidence for a 
first order phase transition when n 1 ' 3 < a em m. In the concluding section we make some 
comments on the validity of the perturbative expansion and the relevance of our findings 
to helium-core white dwarf stars. 

2 Charged condensation at zero temperature 

We start by considering the Lagrangian of a charged scalar field $ of mass m and a gauge 
field at zero temperature. The scalars carry charge g which we take to be some multiple 
of the electric charge g = Ze. We include a background external charge density g' which 
has a charge opposite that of the scalars: 

£ = -~A + \ D ^\ 2 - m2<r$ - • C 1 ) 

The covariant derivative for the scalars is defined as = — igA^. 

This Lagrangian can arise from a gauge invariant theory in which fermions with charge 
g' play the role of the background charge density: 

£ = ~\ F lu + \ D ^\ 2 ~ ™ 2 $*$ + - m F )TP . (2) 

In this work we consider such a theory in the limit that the fermions are non-dynamical 
(i.e., frozen in "by hand" or by some other dynamics). At distance scales that are greater 
than the average separation between the fermions, their spatial distribution can be as- 
sumed to be uniform. Then, the background charge density can be approximated as 

Due to the global U(l) symmetry, there is a conserved scalar current: 

= -*[(£>,,*)** -^(D^)]. (3) 

We can associate a chemical potential /i with the conserved scalar current. To incorporate 
the presence of a conserved current into the Hamiltonian density we make the usual shift 
% — > %' = T~L — ^Jq. For the Lagrangian density the equivalent shift is 

C -> C! = -\F%, + i^o)^ 2 " ™ 2$ * $ " • ( 4 ) 
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For convenience we switch notation and write the scalar field <E> in terms of a modulus and 
a phase <I> = -^o 'e l ° ' . In terms of these variables the Lagrangian density becomes 

£' = -\*% + W^f + \(gA^ - d^a + ^ ) 2 o 2 - \m 2 a 2 - g'A^J^ . (5) 

Written in this form, it is evident that the chemical potential acts as a tachyonic mass for 
the scalar s. 

In these variables the scalar number density is given by 

J S = {gA - d a + n)a 2 . (6) 

Varying the Lagrangian with respect to A^ gives the following equation of motion: 

- d^F^ = g{gA v - d v a + [ib M )a 2 - g' J v . (7) 

The system is electrically neutral when the scalar charge density is equal in magnitude to 
the background charge density: 

g4 = g'Jo ■ (8) 

Varying the Lagrangian with respect to a gives: 

aa=[(gA IM -d fl a + n5f l0 ) 2 -m 2 ]a, (9) 

while varying with respect to a gives the conservation equation for the scalar current: 

d^j; = d» [{gA^ - d^a + iiS^o) a 2 ] = 0. (10) 

Let us take Aj — dja = so that the number density of scalars is constant in time: 
Jq = const. Subject to this constraint, the equation of motion for the scalars ([9]) becomes 

Ua = ^-m 2 a. (11) 



This system has a constant, static solution 

(^) = a/?- ( 12 ) 

V m 

The nonzero expectation value for a indicates that the scalars are in the condensed phase. 
It follows from ([6]) that in the condensed phase {gAo — doa) + /i = m. For an electrically 
neutral system in which the gauge-independent quantity (gAo — doa) is equal to zero, 
condensation occurs when [i = m. Thus at zero temperature we have a neutral system in 
which the charged scalars are condensed into a zero-momentum, macroscopic state. For 
further discussion of the zero temperature condensate see [18] . 
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3 The Thermodynamic Potential 



In order to study the properties of the charged condensate at finite temperature we start 
by computing the grand canonical partition function Z. We use the functional integral 
representation of the partition function as it is most suited to the field-theoretic approach 
adopted above: 



Z = N J [d$p$*][cL4 M ] det (jpj 5(F) exp 



dr I d 3 xC 



o 



(13) 



Here N is an irrelevant normalization constant, /3 = 1/T is the inverse temperature and 
r = it is the imaginary time. We have introduced a gauge fixing condition F in order to 
evaluate the functional integral over the gauge fields. The action is defined by 

S = - I dr [ d 3 xC , (14) 



Jo J 

where the Lagrangian density £ is given by (|4|) and contains the scalar chemical potential 
as well as the external current J^. We also make the replacement Aq = —iAq. 
Let us decompose $ into real and imaginary parts: 

<& = ^(0i+# 2 ). (15) 

The functional integral over the fields is constrained to be periodic so that 0(0, x) = 
(j)(/3,x). Given this constraint, we can Fourier expand the fields as follows: 

fa{x) = v + f; ^e^^*) 0i ;n (p) , (16) 

n=— oo p 

^) = y^EE ei[ " nT +p ' x) ^n(p) , (17) 

n=— oo p 

W = ^EE e i(UnT +P -x) Awn(p) > (18) 
n=— oo p 

where uo n = 2irnT and V is the volume of the system. For the field <pi we have separated 
out a constant part v which is independent of x and r so that 0i ;n =o(p = 0) = 0. Thus v 
represents the thermal average of the field: (cj)i(x)) = v. We do this in anticipation of the 
condensation of the scalars into the n = 0, p = state. The partition function will be at 
a minimum with respect to the free parameter v: 

dlnZ , . 

— = »• < 19 > 

Without loss of generality we can set (02 (x)) = 0. 

One could also separate out an expectation value for the scalar potential: (A§(x)) = clq. 
The partition function will also be at a minimum with respect to a®: 

dkiZ , , 

^— = 0. 20 
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It can be seen from the Lagrangian ([5]) that a nonzero expectation value for the scalar 
potential acts as a shift in the chemical potential, except in its coupling to the external 
current Jo- Thus equation ([20]) is equivalent to 

1 d\nZ . , . 

S w^-- 9J ° =0 - (21) 

The first term on the right hand side is the scalar charge density: 

1 81,12 «■ P2) 



PV dfx 

Thus minimization of lni? with respect to ao gives the condition of charge neutrality: 

g(J s )=g'Jo- 

Equation (I2ip depends only on the total effective chemical potential, i.e., on the sum 
/i e ff = M + gO'O and not on ao alone. Thus minimization does not fix the value of ao 
independently of fi. Instead, equation (|22j) can be used to determine ^ e ff> given a fixed 
charge density. One could consider a configuration with a nonzero ao acting as an external 
chemical potential for the scalars, corresponding to some uncompensated charge on a 
surface at infinity. Such a setup was studied in earlier works [HDS] in the zero-temperature 
limit. For simplicity, we will assume a configuration with no surface charge and take ao = 
in what follows 

Terms in the Lagrangian that are linear in the excitations of the fields above their 
background values will contribute to the action an amount proportional to 0i,2-n=o(p = 0) 
or ^4o;n=o(p = 0) for the scalar and gauge fields respectively. Thus given the above defini- 
tions, these linear terms in the Lagrangian will not contribute to the partition function. 
Terms that are linear in Aj- n (p) do not appear in the Lagrangian as long as the external 
current Jj is zero. 

In order to compute the functional integral we must first choose a gauge F. The most 
natural choice is the unitary gauge in which the phase of the scalar field is set to zero: 
F = a = 0. In this gauge the physical content of the theory is apparent and no unphysical 
degrees of freedom need to be introduced. However, the unitary gauge is known to give 
incorrect results at the one-loop level. This is attributed to the fact that Lagrangian in the 
unitary gauge does not correspond to a renormalizable theory (HJ [9]. We choose instead a 
renormalizable gauge. A family of covariant gauges is given by 

F = + gvfc - f(x) = , (23) 

for some arbitrary function f(x). With this choice the partition function (113p becomes 

Z = N J [#x] [dfa] [dAp] det (□ + <? V) 5(8^ + gvfc - f(x)) e~ s . (24) 

The delta function in the partition function can be incorporated into the Lagrangian via 

£ -+£ -^-(8^ + gvfc) 2 . (25) 



lr To consider a system with a nonzero surface charge, one need only make the replacement jj, —¥ /j, e s in 
the following expressions. The results are otherwise unaffected. 
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In what follows we will keep the gauge parameter £ general in order to check that our 
results are independent of the gauge fixing condition. 



To compute the partition function we proceed via the mean field approximation. We 
take the field expansions (fTB]) . (fT7|) and (fT8|) and substitute them into the Lagrangian 
(125|) . We expand the Lagrangian to second order in the fields, which we assume to be 
small fluctuations above the mean field values. We neglect the terms that are linear in 
the fields, as per the discussion above. For convenience we split the gauge boson into 
transverse and longitudinal components, the transverse components being given by 



4 r ' 



tr 



A, 



V 2 



(d k A k 



(26) 



Upon integrating, the action (to second order) can be written as a sum of the tree-level 
component, the contribution coming from the transverse photons and a contribution com- 
ing from the scalars and remaining photon degrees of freedom: 



S2 — So + Str + S s . 

The tree-level action is given by 

S = \l3V{m 2 - fi 2 )v 2 



(27) 



(28) 



The contribution to the action coming from the transverse components of the gauge boson 
is 



Str = ^EEE 4[- n (-p)K + P 2 + sV)< n (p) , 
n p i=l,2 



(29) 
(30) 



where the sum over % is the sum over both transverse degrees of freedom. These components 
acquire a mass when the scalar is condensed, i.e., when v 7^ 0. 

The remaining scalar and photon degrees of freedom mix with each other. Their 
contribution to the action is given by 



>h-n{ 



-p), 2; _„(-p), Ao--n(-p), ^ ; _„(-p)) D 



n p 



/ 01;n(p) \ 

Ao;n(p) 
\ 4;»(P) ) 



(31) 



where D = 
( + p 2 + m 2 - n 2 

+ 2lJLLO n 

—2igv[i 




—2fj,ui n —2igvji \ 

I + p 2 + m 2 - n 2 + \g 2 v 2 +igvu n (l-h +igv\p\ (1 - j) 

igvu n (l-\) \ujl + p 2 + g 2 v 2 -w„|p|(l-|) 

-u n \p\{l-\) ul + \p 2 +g 2 v 2 j 



-igv\p\{\ - -A 



(32) 



7 



Let us make a few comments on this matrix. In the decoupling limit g — > the scalar 
degrees of freedom decouple from the photon degrees of freedom as expected. In this 
limit we recover the theory for the neutral condensate together with decoupled, massless 
photons. When g ^ but the scalar field is not condensed, i.e., when v = 0, the same 
decoupling takes place. In the uncondensed phase the partition function for the charged 
scalars is the same as that for the neutral scalars and decoupled photons, at the one-loop 
level. 

An important element of this matrix is the term that mixes the scalar field (p\ with 
the 0-component of the gauge field: 2igvfMpiAo. In the absence of a chemical potential 
for the scalars (i.e., when // = 0), an appropriate choice of gauge will decouple the gauge 
degrees of freedom from the scalars, even when v is nonzero. When this decoupling 

no longer occurs, because of this term. Due to this mixing, the spectrum of the theory in 
the condensed phase is significantly different from that of the neutral condensate. 

Carrying out the functional integrations, the logarithm of the partition function (|24ft 
can be expressed as a sum of the various components 

In Z = In Zq + In Zt r + In Z s + In Z g h . (33) 

There is also a constant term coming from the overall normalization which is temperature 
independent and which we can neglect. 

The first term comes from the tree-level action 

lnZ = -\f3V{m 2 - ^i 2 )v 2 , (34) 

the second term from the two transverse photon polarizations 

lnZ tr = -i(2)/3 2 Y, E lndet (^n + P 2 + 9 2 v 2 ) , (35) 

n p 

and the third term from the scalar and remaining photon degrees of freedom 

lnZ s = -±/3 2 ^^ ln(£det£>). (36) 

n p 

The last term arises from the ghost determinant in the partition function (|24p . 

In Z gh = (3 2 Y E ln det ^n + P 2 + 9 2 v 2 ) ■ (37) 

n p 

Formally, the ghost term exactly cancels the contribution coming from the transverse 
degrees of freedom. This allows us to simplify our calculations. 
The determinant of D can be written in the following form 

£ det D = (u 2 + u 2 + )(uj 2 n + u 2 _)(u 2 + u 2 ){u 2 n + u 2 ) , (38) 

where the dispersion relations w+, w_, u)\ and 0J2 are functions of p, m, fx, gv, and £. We 
will denote them by w a . Due to their length we will not write the full expressions here. 
Below we will give their exact expressions when evaluated on the solutions to the tree-level 
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equations of motion for the scalar. We will also find expressions for when the tree-level 
scalar equations of motion are not satisfied. 

This form for £ det D allows us to carry out the sum over n in (|36p (see [7j for more 
details). We replace the sum over p by the appropriate integral. The thermodynamic 
potential = —-^y InZ is then: 

n = \{m 2 - ,V +E / + / ■ ^ 

The first term is the tree-level potential. The second term is due to the zero point energies 
of the fields. The third term comes from the thermal excitations of the fields. 

The thermodynamic potential is subject to two constraints. First, it should be at a 
minimum with respect to the free parameter v. 

For the tree-level potential Ot ree = \{m 2 — [i 2 )v 2 we see that the first condition is satisfied 
in two ways: in the condensed phase when v ^ and fi = m or in the uncondensed phase 
when v = and \i 7^ m. We will refer to these solutions as the "on-shell" solutions since 
they satisfy the tree-level equations of motion for the scalar. 

Second, the thermodynamic potential should give the fixed number density of scalars 
when differentiated with respect to /z: 

»-W=-f. (41) 
When applied to the tree-level potential, this constraint gives 

n = ^— = fj,v z . (42) 



In the condensed phase this is satisfied when v = Wn/m, as we found in the previous 
section for the condensate at zero temperature. In the uncondensed phase (v = 0) this 
expression cannot be satisfied. This is unsurprising as we have neglected contributions to 
the number density coming from thermal fluctuations. In the absence of these fluctuations, 
i.e. at zero temperature, the charge of the system must be stored in the condensate. 

The dispersion relations u a simplify greatly when evaluated on the on-shell solutions. 
For the uncondensed phase we set v = in the determinant of D (equation (|32p ). Then 
det-D can be factored according to equation ([38]) . We find: 



oo± = ^p 2 + m 2 ± fi , wi )2 = |p| • (43) 

The first two dispersion relations are for the scalar particle and antiparticle and are the 
same as those for a neutral system at finite chemical potential. The second two are for 
the photon degrees of freedom. In the uncondensed phase, the spectrum of the theory 
and thus the one-loop thermodynamic potential is the same as for neutral scalars and 
decoupled, massless photons as mentioned above. 
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To find the dispersion relations in the condensed phase we set \x = m in the determinant 
of D. This gives 

u;| = P 2 + \g 2 v 2 + 2m 2 ± y'Wp 2 + (2m 2 - ±# 2 ^ 2 ) 2 , (44) 



< 2 = P 2 + 3 V • (45) 

These dispersion relations agree with those found for the charged condensate at zero tem- 
perature in [I], with the photon mass m 7 = g^n/m replaced by the more general gv. In 
the condensed phase the first two dispersion relations (I44D no longer have the simple inter- 
pretation of corresponding to the scalar particle and antiparticle. In the condensate the 
gauge symmetry is broken. The gauge boson becomes massive by eating one of the scalar 
degrees of freedom. Thus we can think of these two dispersion relations as corresponding 
to the longitudinal component of the photon and the remaining scalar degree of freedom, 
though in reality they correspond to a linear combination of the scalar and gauge fields. 
Their masses are found by taking p = 0: w±(p = 0) = 2m, gv. Both modes are massive. 
Thus unlike the neutral condensate which contains a massless particle in the condensed 
phase, the charged condensate has a mass gap. 
For m^> gv the relations (I44D simplify to 



~ \/p 2 + m 2 + m , (46) 



and 



,2/V»2 „2„,2> 



2 2 2 , P \P 9 v ) 2 ^ n i »m\ 

ujI ~ g z v z H — -= , p z < , (47) 

w_ ~ \/p 2 + itt- 2 — m, p 2 ^ 2gvm . (48) 



The solutions to equation (j4*0|) calculated from the tree-level potential will be modified 
by finite temperature effects. In particular, we will find that at non-zero temperature, 
\i = m no longer holds identically in the condensed phase. Moreover, in order to solve 
the second constraint equation (|4ip at finite temperature, we need to know the dispersion 
relations oj a as a function of fj,, away from \x = m. Thus to understand the full thermo- 
dynamic potential and its constraints, we must generalize the u a given above to the case 
that the tree- level solutions to the scalar equations of motion are not satisfied, when v ^ 
and fj, 7^ m. 

To compute the thermodynamic potential we have used the standard background field 
method. For gauge theories this method gives a unique result only when the background 
field is a solution to the tree-level equations of motion. If this condition is not satisfied, 
the background field method can give gauge dependent results for physical quantities. 
This effect can be seen here by considering the matrix D given in (I32p . Calculating the 
determinant of D, one finds terms that depend on the gauge fixing parameter in the form 
£(m 2 — ^ 2 )g 2 v 2 {. . .). These terms vanish on-shell, i.e., when either v = or fi = m. Thus 
the dispersion relations found above are independent of £ as we would expect. However, 
when d / and fi ^ m the determinant of D and thus the more general dispersion 
relations become dependent on the gauge fixing condition. 
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In order to determine the general, gauge-independent dispersion relations from det D 
we would first have to find the unique off-shell potential. For our purposes, however, it is 
possible to determine the more general uj a , given our knowledge of their behavior when 
v = or [i = m, in combination with constraints coming from a Ward identity. We will 
do this now. 

The photon degrees of freedom do not contribute to the number density of the scalars. 
Thus the generalized dispersion relations oj\^ must be independent of [i. It follows that 
they are the same when /i ^ m as when [i = m: 

ul 2 = p 2 + g 2 v 2 . (49) 

To find the dispersion relations for the scalar degrees of freedom, we start by parametrizing 
(A by 



U) 



l = A±VB. (50) 



By requiring that the general dispersion relations reduce to those found above when v 
or when fj, = m we find 



l 2,2 2 1 2 2 1 / \ "' /' 2 2 

A = p+m+fi+-gv +a{m,gv,LL) = — g v 

1 



B = Ap 2 fi 2 + Am 2 f^ -2m 2 g 2 v 2 + -g 4 v^ + b(m,gv,fj,)(m 2 - fj, 2 )g 2 v 2 . (51) 

The functions a and b are not fixed by this requirement. They can in general represent a 
series expansion in m 2 — fi 2 and g 2 v 2 with undetermined coefficients. We take a and b to 
be independent of |p|, in order to maintain the appropriate high-momentum behavior of 
the theory. 

For the neutral condensate, \i = m identically in the condensed phase, even at finite 
temperature. For the charged condensate this is not necessarily the case. This is due to 
the fact that for the charged condensate the dispersion relations depend on v. However, 
any deviation of \x away from m should be suppressed by a em and vanish in the limit that 
g —7- 0. Thus for the low density systems that are the primary interest of this paper, the a 
and b terms are subdominant as they suppressed compared to the leading order 

terms in A and B. For high density systems they can be relevant, when gv 3> m. 

A Ward identity can be used to fix a and b. The sum over the zero point energies in the 
thermodynamic potential (j39H is UV divergent and thus must be renormalized. Because 
of the Ward identity associated with the conserved scalar current, one can show that the 
conserved current is not subject to infinite renormalization, even in cases of spontaneous 
symmetry breaking (see, e.g., [19]). It follows that the divergences in the potential should 
be independent of \i. Expanding the dispersion relations (I50p for large p we find: 

(52) 



f d 3 p [ f 2 (m,gv,p) U(m,gv,/j,) f 1 \ 

J (^F (w++a; - ) " |p| r + — ^ — + — ^ — +0 {^) 



where f 2 and depend on a and b. The three terms on the r.h.s. represent the quartic, 
quadratic and logarithmic divergences of the scalar sector. These divergences can be made 
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independent of /i for appropriate choices of a and b. This uniquely fixes a 
The generalized dispersion relations are thus: 



and b = -2. 



4 



p 2 + m 2 + n 2 + -g 2 v 2 ± 



4p 2 // 2 + Am 2 jx 2 + 2g 2 v 2 fi 2 — 4g 2 v 2 m 2 + ~^9 4 v 4 • 



(53) 



Let us now treat the remaining, ^-independent divergences in the thermodynamic 
potential. To regularize the divergent terms we introduce an ultraviolet cutoff A c and take 
|p| < A c . After integration, and in the limit of large A c , the divergent contributions to 
the thermodynamic potential coming from both the scalar and gauge degrees of freedom 
are 



d 3 p 

W) 



3 2 



A 4 



Air 2 



+ 



m 2 k 2 c 



+ 



8vr 2 
3A 2 3m 2 



m 

W^ 2 



In 



167T 2 
1 

32^ 



—A— 



2A r 



m 



2 2 

9 v 



In 



2A f 



m 



+ 



1 



16vr 2 



In 



2A C 
gv 



g A v 4 + finite terms . 



(54) 



The first term renormalizes the vacuum energy density. This term is independent of 
temperature and can be trivially subtracted. The second term renormalizes the mass of 
the scalar and the third term renormalizes the coupling of the cj) 4 interaction, had we 
included one at the start. We can add to the Lagrangian counterterms of the form 

5m 2 ®*& 



-c.t. 



<5A 



(55) 



to absorb these divergences. 

To renormalize the zero-temperature thermodynamic potential to leading order in a em , 
we can evaluate the contributions from the zero-point energy on the tree-level solution 
= m. We also, for the moment, ignore the divergences that multiply v 4 as they are also 
higher order in a em . We choose our renormalization conditions so that at zero temperature 
the thermodynamic potential is finite and is at a minimum when with respect to v when 
li = m. This fixes 5 A and 5m 2 so that 

+ ^z.p. + ^c.t. = \{m 2 - fi 2 )v 2 + 0(a 2 cm ) . (56) 

The presence of a v 4 term in the potential will shift the zero-temperature values of v and 
/i. Moreover, the dependence of the divergent term that multiplies v 4 on the logarithm of 
v means that the coupling of the v interaction will run mildly with temperature. However 
this term is suppressed by a em compared to both the other zero temperature terms and the 
finite temperature terms. We will neglect its contribution to the potential when studying 
finite temperature effects in what follows. 

4 Phase Transitions 



In section 3 we found that the minimization of the tree-level potential with respect to v has 
two solutions: one when v = and one for arbitrary v when \x = m. We also observed that 
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at zero temperature the uncondensed solution v = cannot satisfy the constraint equation 
()41|) . At zero temperature all of the scalars reside in the condensate. More generally, there 
exists a critical temperature T cl below which some fraction of the conserved scalar charge 
must be in the condensed phase. Likewise, there exists a critical temperature T C 2 above 
which v = is the only solution to both equations ([40 p and (|4ip . Between these two 
temperatures the scalar field will undergo a phase transition into the condensed state. In 
this section we will determine T c \ and T C 2 by applying the constraints (|40p and (|4ip to 
the full thermodynamic potential. We will also discuss features of the phase transition. 

Let us consider the nonrelativistic limit n 1 / 3 <C m. In this case w+ ~ 2m. Since 
gv <C m in this limit, the cj+ term in the thermodynamic potential is exponentially 
suppressed compared the uj- and uj\ 2 terms. Thus we can neglect its contribution: 



i(,„ 2 - ,V + r / ^ m (1 - .-*-) + 2r/ *t m (1 - e -«v^W) . (57) 



The thermodynamic potential should be at a minimum with respect to u: d£l/dv = 0. 
Since t> appears in the dispersion relations only in the form g 2 v 2 we can write this condition 

as 



(m 2 — fj, 2 )v + 2g 2 v 



d 3 p I duj- 11 1 



(2vr) 3 \d(g 2 v 2 ) eP"- - 1 ^ p 2 + g 2 v 2 6 /?Vp^W _ 1 



(58) 



This equation has two solutions. The first solution is the same as for the tree-level poten- 
tial, when v = 0. The other solution is given by 

2 _ 2 2 f d 3 p ( doo- 1 1 1 \ 

M " m +25 J (2^ - 1 + TJl + 5 2^2 e/?v /?W _ J • (59) 

At finite temperature, /x is shifted away from its tree-level value m by an amount propor- 
tional to a em . It is straightforward to see that w_ as given by ([53]) increases monotonically 
as a function of g 2 v 2 , i.e., dco-/d(g 2 v 2 ) > for arbitrary v. Thus at finite temperature, 
the second solution to d£l/dv = requires jx> m. 

For the neutral condensate, the integral in the thermodynamic potential corresponding 
to the ui- mode is convergent only when /x < m. For the charged condensate however, 
when v ^ 0, solutions with fi > m are possible because of the contribution gv makes to 
the mass of this mode. Still, when v = we must have /x < m, as in the neutral case. 

It follows that v = is not a solution to the above equation (159p . The absence of 
a second solution at v = implies that v must change discontinuously when going from 
the uncondensed phase to the condensed phase. This is indicative of a first order phase 
transition. We will see that this discontinuity is more pronounced in the low density 
regime. 



Let us consider the v = solution. The general dispersion relations given by ([53 
reduce to those found above in ([43jl . In particular 



vV + m 2 - 11 . (60) 
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Figure 1: The scalar expectation value v as a function of temperature for various values 
of n 1 / 3 /a em m when fi = m and g = 2e. 

Differentiating the potential with respect to \x gives the usual result for number density: 

dn = r d 3 P 1 

dfjL J (2vr) 3 e/3 (V^W-^) _ i ' 

All of the particles are in thermally excited states. Fixing n gives an implicit expression 
for the chemical potential as a function of temperature. The critical temperature T c \ is 
defined to be the minimum temperature at which all the particles are still in excited states. 
It is found by taking fj, = m and solving for T. The result is the same as for the critical 
temperature of the neutral condensate: 

T — for 1,1/3 « m - (62) 

The temperature T c \ puts a lower bound on the condensation temperature. Below this 
temperature no v = solution exists. However, unlike in the neutral case, it may be 
possible for the charged scalars to condense before T c \ is reached, as we will now discuss. 

When j; / the number density is given by 

2 f d 3 p &j_ 1 

n = ^ v -J(2^^T^— V (63) 

where cj_ is given by (|53p . The first term on the right hand side represents the fraction 
of particles in the condensate while the second term gives the particles that remain in 
thermally excited states. This equation combined with (|59p gives an implicit expression 
for the value of the condensate v as a function of temperature for fixed n. As noted above, 
the finite temperature correction to /u 2 — m 2 is suppressed by a em . As a lowest order 
approximation we can evaluate the above expression (|63[) on the tree-level solution fj, = m 
to find v(T). 

In Figure 1 we plot v as a function of temperature for various values of n 1//3 /a em m 
using expression (I63p and taking fj, = m. In the relativistic regime we also include the 
contribution to number density coming from the w+ mode. We take the charge of the 
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scalars to be g = 2e with a system of helium nuclei in mind. The values of v at T = and 
of the critical temperature T c \ are of course different for different values of m and n, but 
we rescale the results for easy comparison. 

In Figure la we plot v(T) in the nonrelativistic regime when n 1 / 3 < m. The orange 
line corresponds to n 1 / 3 = 1/10 a cm m, the red line to n 1 / 3 = a cm m, and the pink line to 
n 1 / 3 = 10 a cm m. For comparison we have also plotted v(T) for the neutral condensate in 
the nonrelativistic limit. This is the black dashed line. In Figure lb we plot v{T) in the 
relativistic regime when n 1 / 3 > m. The purple line corresponds to n 1 / 3 = 10 2 a cm m, the 
dark blue line to n 1 / 3 = 10 3 a cm m and the light blue line to n 1 / 3 = 10 4 a em m. The black 
dotted line is v(T) for the neutral condensate in the relativistic limit. All of the curves 
share the solution v = for T > T c \ . 

The shape of the contours for n 1 / 3 < a cm m is characteristic of a first order phase 
transition. At sufficiently high temperatures the only solution is v = 0. But as the 
temperature drops, at some temperature T C 2 > T c i, the v(T) curve ceases to be single- 
valued. Both the uncondensed solution and the condensed solution exist simultaneously. 
Moreover, in between and T c \ the value of v must change discontinuously from zero 
as the system cools. At temperatures below T c \ only the condensate solution exists. 

The temperature T C 2 is the maximum temperature at which the condensate solution 
exists. It corresponds to the point on the above plots where 



dv 
df 



= oo . (64) 

T=T c2 



We see from Figure 1 that lowering the value of n 1//3 /a em m appears to increase T C 2 relative 
to T c i. Lower densities also appear to correspond to larger discontinuities in v. 

For the charged condensate, as n 1//3 increases relative to a cm m the function v(T) 
appears to asymptotically approach a single curve that is similar to the curve for the 
neutral condensate. The phase transition appears to be nearly very second order with 
v increasing more-or-less continuously as T decreases. Even though the phase transition 
may technically be first order the jump in v becomes negligibly small. 

The above contours were drawn using the approximation /j, = m in the condensed 
phase. The actual condensate solutions are slightly to the left of these contours, when 
^ = m + 0(a cm T). 

The existence of a first order phase transition at low densities can be understood 
in the following way. In the low density regime the mass of the bosonic excitation gv 
becomes comparable to or greater than the transition temperature. Thus the thermal 
contributions of this mode to the number density become exponentially suppressed by 
a factor of exp(— gv/T) in the condensed phase. In order to maintain the fixed charge 
density, v must jump discontinuously so that charge that had been stored in thermal 
excitations is stored instead in the condensate. As the density is lowered, T c continues to 
decrease relative to gv and the effect grows more pronounced. 



as 



Let us now find T C 2 . In the nonrelativistic limit we can approximate the number density 

2 , f d3 P P 2 1 ra^ 

n~m,V + / rr a > (OOj 

J (2vr) 3 2mw_ e^- - 1 
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Figure 2: Critical temperatures T C 2/T C \ as a function of a em m/n 1//3 when g = 2e. 



^p^/Am 2 + g 2 v 2 . Here we have taken /i 



where ex- 
changing to dimensionless variables 



m in the condensed phase. 



p- 



2mT 



and 



the number density becomes 



- . ( mT \ 
\2n J 



3/2 



gv_ 
T 



dx x 3 / 2 



7T Jo ^x 2 + y 2 e V x2 +v 2 



(66) 



(67) 



Let us assume that at the critical temperature T C 2 the expectation value of the con- 
densate as given by the above expression is roughly the same order of magnitude as it is 
at zero temperature: v ~ \Jnjra. Let us also assume that T C 2 is roughly the same order of 
magnitude as T c \ as given by (|62p . Then when re 1 / 3 <C a cm m we have T C 2 <C gv. In terms 
of our dimensionless variables this is y S> 1. In this regime we can do a low temperature 
expansion of the above integral |21| . The first two terms of the low temperature expansion 
give: 



ii 



mv + 



2 5 / 4 T[5/4] fmT\ 



3/2 



7T 



gv\V 4 



Li. 



L 5/4 



(-'" /T )-|(?)" 3/1 



Li 



L 9/4 



(e"^/ T ) 



where Li is the polylogarithm. 

T C 2 corresponds to the maximum temperature at which this expression can be satisfied 
for a fixed n. Equivalently, it corresponds to the critical point of the function v(T). So to 
find T C 2 we differentiate both sides of the above expression with respect to v and we set 
dT/dv = 0. We find that at the critical temperature, T C 2 and v are related by 



V C 2 




(69) 



Substituting this expression for v into (I68j) gives an implicit expression for T C 2 in terms of 
n and m. 
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In figure 2 we plot the ratio of T C 2 to T c \ as a function of aemm/n 1 / 3 , from aemm/n 1 / 3 = 
1 to 500. Again we set g = 2e. As n decreases the ratio T C 2/T C \ grows. Since T c \ coincides 
with the critical temperature for the neutral condensate, this implies that for the charged 
condensate at low densities, condensation can occur at temperatures several times higher 
than for the neutral condensate. In high density regimes where a era m/n 1 ^ < 1, the above 
approximations break down. We expect that in this limit the Coulomb energy of the scalars 
should be small compared to their kinetic energy and thus T C 2/T C \ should approach 1. 

5 Conclusions and Discussion 

We have calculated the one-loop grand canonical partition function for scalar electrody- 
namics at finite temperature and chemical potential. We considered an electrically neutral 
system in which the charged scalars reside in a background density of the opposite charge. 
Using the background field method we found a result that was gauge condition dependent 
when the background fields did not satisfy the tree-level equations of motion. To address 
this, we required consistency with a Ward identity in order to determine the off-shell 
thermodynamic potential. 

Using this potential we could determine, to lowest order in perturbation theory, the 
thermal expectation value of the charged scalars v as a function of temperature. Plotting 
v(T) for various values of the scalar mass m and number density n, we found evidence of 
a first order phase transition in the low density regime, when n 1 / 3 < a em m. In between 
two critical temperatures T c \ and T C 2, v jumps discontinuously from zero. The lower 
temperature T c \ coincides with the critical temperature for a neutral condensate, while 
T C 2 can be several times higher in the low density regime. Thus the phase transition for 
the charged condensate can occur at higher temperatures than for the neutral condensate. 
The strength of the first order phase transition appears to increase with decreasing density. 

The critical temperatures T c \ and T C 2 give only lower and upper bounds on the tran- 
sition temperature. Finding the exact transition temperature as well as other physically 
interesting quantities, such as the latent heat involved in the phase transition or the 
timescale for the transition to occur, requires a more thorough understanding of the ener- 
getics of these solutions and their stability or metastability. Investigation of these questions 
would entail using the constraint equation n = —d£l/d[x to "integrate out" /x so that the 
free energy density F(fi,v,T) = fi(^,u,T) + fj,n — > F(n,v,T). This is analogous to our 
procedure for the zero-temperature potential in equation (|lip . 

The above conclusions were based on a one-loop calculation. It is reasonable to consider 
the effects of higher order corrections in perturbation theory. In thermal field theory at 
high temperatures, the perturbative expansion can break down, even for weakly coupled 
theories. This can occur when the temperature of the system is greater than the masses 
of the particles running in the loops. For theories with massless particles, this breakdown 
is often signaled by the appearance of infrared divergences at higher loops. 

For the system described here, because the gauge symmetry is spontaneously broken 
in the condensed phase, the gauge bosons are massive. Thus one can compute to any 
order in perturbation theory without encountering infrared divergences. Moreover, in the 
condensate, as long as the temperature is sufficiently below the gauge boson mass gv and 
the mass of the remaining scalar degree of freedom 2m, the loop expansion parameter 
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v(T) 




Figure 3: Region of validity of the perturbative expansion 

should remain small. In the uncondensed phase, the massless photons interact with each 
other only through the massive scalars, the lightest of these scalars having mass m — Li. 
Thus as long as m — fi S> T our considerations in the uncondensed phase should be valid. 

For a second order phase transition or for a weakly first order phase transition these 
criteria will be violated near the phase transition point, where gv — > and fx — )■ m. Thus 
perturbation theory can not be considered reliable in the vicinity the phase transition. 
However, for a strongly first order phase transition, which is precisely the interest of this 
work, this region is avoided. In the low density regime the condensate solution doesn't 
pass near gv = 0; the temperature of the system is always small compared to the masses 
of the excitations in the condensate. Thus higher order corrections to the results found 
above should be small. 

To illustrate this we plot v(T) in figure 3 for the values of m and n that are relevant 
for helium-core white dwarf stars (m ~ 4 x 10 4 n 1//3 ). The orange region corresponds to 
where T > gv and thus perturbation theory can be unreliable. The condensate solution 
which is the upper solid blue curve is well outside this region. We also mark the point 
where, when v = 0, T = m — fi. To the left of this point, i.e., when T > m— fx, calculations 
of, say, the free energy density in the uncondensed phase using the expressions obtained 
above should not to be trusted. This may place some limitations on calculating the exact 
transition temperature, if the phase transition occurs in this region. However, it is possible 
that methods such as resummation can be applied in this region to obtain more reliable 
results. 

Let us end with some comments on physical applications. In helium-core white dwarf 
stars the number density of the helium nuclei is such that n 1 / 3 < a em m. Thus based on the 
arguments given above, we might expect the condensation of these nuclei to be a first order 
phase transition. We may also expect the condensation temperature for the helium nuclei 
to be somewhat greater than that for neutral bosons. The latent heat associated with the 
phase transition could potentially delay the cooling of the star. For a carbon-core white 
dwarf star of mass ~ Mq which crystallizes, the latent heat associated with crystallization 
can increase the classical cooling time by a factor of ~ 1.6 [22] . However, one of the main 
features of a condensed-core white dwarf is that it cools much more rapidly than white 
dwarfs with uncondensed or crystallized cores [5]. It is possible that any contribution to 
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the cooling history coming from the latent heat will be overwhelmed by the rapid cooling 
in the condensed phase. 

In this work we treated the background charge density as static and classical. For the 
system of nuclei and electrons in the cores of white dwarfs, this is not a good approxima- 
tion. The electrons, i.e. the background charge density, form a degenerate fermi gas with 
gapless excitations. Due to these excitations, the electrons could play a significant role in 
the phase transition. 
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